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One of the main limitations of utilizing optimal wavefront shaping in imaging and authen¬ 
tication applications is the slow speed of the optimization algorithms currently being used. 
To address this problem we develop a micro-genetic optimization algorithm (/iGA) for opti¬ 
mal wavefront shaping. We test the abilities of the /xGA and make comparisons to previous 
algorithms (iterative and simple-genetic) by using each algorithm to optimize transmission 
through an opaque medium. From our experiments we find that the pGA is faster than 
both the iterative and simple-genetic algorithms and that both genetic algorithms are more 
resistant to noise and sample decoherence than the iterative algorithm. 


I. INTRODUCTION 

The ability to use wavefront shaping to con¬ 
trol the optical properties of opaque media was 
first predicted by Freund in 1990 [1] and demon¬ 
strated experimentally in 2007 by Vellekoop and 
Mosk; who used a liquid crystal on silicon spatial 
light modulator (LCOS-SLM) to shape the wave- 
front of a laser such that the beam was focused 
through the material [sj]. Since then the tech¬ 
nique of using optimal wavefront shaping in con¬ 
junction with a scattering material has been used 
for polarization control Q, 4], spectral control 
it 5-^], enhanced fluorescence microscopy 


niques have been developed to aid in wavefront 
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photoacoustic wavefront 


One of the key features of all of these appli¬ 
cations is the use of an optimization algorithm 
to search for the optimal wavefront. The choice 
of which algorithm to use is crucial and deter¬ 
mines the speed and efficiency of optimization 
as well as the algorithm’s resistance to noise and 
decoherence effects [26H29I]. In this paper we will 


lOf] , perfect focu sing {s, 11], compression of ul¬ 
trashort pulses [12, 13], spatio-spectral control of 


ranced astronom- 
19! ]. Also new tech- 


focus on the software (algorithm) used in opti¬ 
mization. However, we note here that the hard¬ 
ware used also plays a major role in determining 
the speed of optimization, as wavefront optimiza¬ 
tion is fundamentally limited by the refresh rate 
of both the SLM and feedback detector. Thus far 
there have been three main classes of algorithms 
used in the literature: iterative 12, 26], parti¬ 


tioning 


26], and simple-genetic 4 , \2m\- While 
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each method is found to have different benefits 
and limitations, one of the common limitations 
is slow optimization speed. In order to obtain 
faster optimization speeds we develop a micro- 
genetic algorithm (^rGA) for wavefront optimiza¬ 
tion as /iGAs are known, from other applications 
fl 31], to be faster than both the iterative and 
simple-genetic algorithms. 

In this study we first describe, in detail, the 
operation of the iterative, simple-genetic, and 
micro-genetic algorithms. We then discuss the 
basic differences between the algorithms and how 
these differences lead to the micro-genetic algo¬ 
rithm being the fastest algorithm. Finally we 
perform an experimental comparison of the three 
algorithms using an optimal transmission experi¬ 
ment [29( to test the algorithms’ speed, efficiency, 


is systematically evaluated in order to find the 
optimal solution [32]. While this technique is 


guaranteed to work in the absence of noise, it is 
extremely time consuming - becoming unfeasi¬ 
ble as the number of dimensions in the solution 
space increases - and ineffective as the signal-to- 
noise ratio decreases |27‘,, 28]. 

A different approach to optimization prob¬ 
lems, which draws inspiration from nature, are 
so called evolutionary algorithms (EAs) [33 351]. 


EAs are stochastic algorithms in which solutions 
are randomly generated, evaluated, and then 
modified until the best solution is found. One of 
the most well known classes of EAs are genetic 
algorithms (GAs) and they function based on the 


same principles as natural selection 


33j, 


34, 


36- 


and resistance to noise and decoherence. 

II. THEORY 

A. Background 

Focusing light through opaque media is in¬ 
herently an optimization problem in which we 
seek to find a specific solution (phase pattern) 
which optimizes the evaluation of a function 
(light focused in a given area). Optimization 
problems have been extensively studied and nu¬ 
merous techniques have been developed to ap¬ 


proach optimization 


32l] . One of the simplest 


approaches to optimization problems is the brute 
force method in which every possible solution 


381 ]. In GAs a population of solutions is ran¬ 
domly generated and then evaluated to deter¬ 
mine how well each member solves the optimiza¬ 
tion problem. The functional response of eval¬ 
uating the i th population member is called the 
member’s fitness, /,;. Once the members of a 
generation are ranked by their fitness scores, a 
new generation is produced through a process 
called crossover. During crossover the algorithm 
chooses two different “parent” solutions using a 
fitness based selection method and then the two 
parents are combined such that the new “child” 
solution contains half the values of one parent, 
and half the values of the other parent, similar 
to biological reproduction. This process repeats 
until a new generation of solutions is produced 
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and then the algorithm restarts with the new 
generation. 

The standard implementation of GAs has al¬ 
ready been applied to light transmission opti¬ 
mization [4, 27] with impressive results, espe¬ 
cially in regards to the algorithm’s resistance to 
noise j^T:]. While the standard-genetic algorithm 
is an improvement over the brute force itera¬ 
tive method, it still requires a large number of 
function evaluations to work, which is time con¬ 
suming. An alternative, and faster, GA imple¬ 
mentation is the micro-genetic algorithm (/xGA) 
{sol . 31]. The micro-genetic algorithm primarily 
differs from the standard-genetic algorithm in its 
population size, required use of elitism (carrying 
over the most fit population from a generation to 
the next), crossover technique, and use of popu¬ 
lation resets instead of mutation. In the follow¬ 
ing sections we will describe the SGA and /xGA 
(as well as the iterative algorithm) in detail and 
expand on the differences between the two ge¬ 
netic algorithms. 


B. Iterative Algorithm 


The simplest optimization algorithm for con¬ 
trolling transmission is the iterative algorithm 
(IA) [2, 26]. To begin the IA bins the SLM pix¬ 


els into N bins and then, starting with the first 
bin, the IA changes the bin’s phase value in M 
steps with a step size of 2ir/M until the phase 
value giving the largest target intensity is found. 


This optimal value is then set for the bin and 
the procedure continues through all the bins un¬ 
til a phase front giving peak intensity is found. 
This optimization scheme is simple to implement 
(requiring minimal coding) and accurately inves¬ 
tigates the solution space. On the other hand, 
the algorithm requires M x N function evalua¬ 
tions, which quickly becomes impractical as the 
number of bins increases, and the IA’s efficiency 
is severly limited by the persistence time of the 
signal-to-noise ratio of the system 


sample and 
3, 3, 29]. 


C. Simple Genetic Algorithm 

The first GA we use is a simple-genetic algo¬ 
rithm, which is similar to those previously used 
by Conkey |3] and Guan P] . The algorithm 
proceeds as follows: 

1. Generate 30 binned random phase masks 
using a random number generator. 

2. Evaluate the fitness of each phase mask. 

3. Pass the top fifteen phase masks with 
highest fitness to the next generation and 
discard the bottom fifteen. 


4. Randomly choose two members of the 
fittest fifteen ( A and B where A ^ B) with 
the probability of being chosen given by: 

fi 


Pi = 


w 

ic /j 
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w 

where J2 fj is the net fitness of the popu- 

j 

lation. 

5. Combine A and B into a new phase pat¬ 
tern C, where each bin in C has a 50% 
chance of being the value from A and a 
50% chance of being from B. 

6. During crossover, at each bin in A and B, 
the two parents bin values are compared 
and if they are the same, a running counter 
is incremented. After each pair of bins has 
been compared this counter contains the 
number of bins that have the same value 
between the two parent images. This value 
is then divided by the total number of bins 
giving the similarity score. If the simi¬ 
larity scores of all parents in a generation 
are greater than 0.97 the bottom 15 phase 
masks are dropped and randomly regener¬ 
ated, while the top 15 are kept. 

7. After crossover and the similarity test, the 
bins in C are mutated at a rate of 0.005, 
which means that a given bin has a 0.5% 
change of being assigned a new random 
value. 

8. Steps 3-5 are repeated until the last fif¬ 
teen members of the new population are 
created. 

9. The fitness of the new population is mea¬ 
sured and steps 3-6 are repeated until a 


stop criteria is reached (set fitness or set 
number of generations). 

D. Microgenetic Algorithm 

The second class of genetic algorithms we 
consider are micro-genetic algorithms (/rGAs). 
H GAs have a similar structure to SGAs, but are 
designed to work with smaller population sizes 
and therefore require fewer function evaluations 
than SGAs. For instance, we typically use a pop¬ 
ulation of 30 phase masks for the SGA, and 5 
phase masks for the /zGA. 

Our ^GA is structured as follows: 

1. Generate five binned random phase masks 
using a random number generator. 

2. Evaluate the fitness of each phase mask 
and rank them 1-5, with 1 having the high¬ 
est fitness and 5 having the lowest fitness. 

3. Discard the two phase masks with the low¬ 
est fitness scores and pass the highest- 
scoring phase mask to the next generation. 

4. Then perform crossover on the top three 
ranked phase masks using the following 
combinations: rank 1 with rank 2, rank 
1 with rank 3, and then rank 2 with rank 
3 is performed twice. 

5. During crossover, at each bin in the par¬ 
ent phase masks, the bin values of both 
parents are compared and if they are the 
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same, a running counter is incremented. 
After each pair of bins has been compared 
this counter will contain the number of 
bins that have the exact same value be¬ 
tween the two parent images. This value 
is then divided by the total number of bins 
giving the similarity score. If the similar¬ 
ity scores of all parents in a generation are 
greater than 0.97 the most fit phase mask 
is kept and the bottom four phase masks 
are randomly regenerated. 

6. Evaluate the fitness of the new generation. 

7. Repeat steps 2 through 6. This procedure 
loops until a stop criteria is reached (set 
fitness or set number of generations) 

E. Differences between the SGA and /iGA. 

While the structures of the SGA and the 
/iGA are similar - both using randomly gener¬ 
ated populations and crossover to explore the 
solution space - there are four main differences 
between the two algorithms: 

1. The /iGA works with a smaller population 
size (five for our algorithm) than a SGA 
(our SGA uses thirty) 


portionate selection 
versa 



39], stochastic uni¬ 


sampling 41], tournament selec- 


40( and reward-based selection 
43j ] . Our SGA uses fitness proportion¬ 


ate selection. 


3. A /iGA requires the use of elitism, while 
an SGA does not require it. 

4. SGA’s use mutation in order to avoid 
local maxima in the fitness landscape, 
while /iGA’s do not. /iGAs rely entirely 
on restarting every time the population 
reaches a local maxima (where the simi¬ 
larity score is greater than 0.97). This is 
step 5 in the /iGA. 


The smaller population size (and lack of mu¬ 
tation) allows the /iGA to optimize more quickly 
[dol . 31 ] as each generation has fewer function 
evaluations than the SGA and the /iGA skips 
the mutation step. However, since the popula¬ 
tion size is smaller, the /iGA requires the use 
of elitism to operate and continue to approach 
an optimal solution 31]. Additionally, the 
small population size and lack of mutation leads 
to the /iGA often drifting toward local maxima, 
which requires the algorithm to restart in order 
to escape a localized solution. 


2. When choosing parents to cross over a 
.ses a tournament style selection 
, while an SGA can use a variety 
of different methods, such as fitness pro- 



III. EXPERIMENTAL METHOD 

For testing the different algorithms we use 
two different types of opaque media: comrner- 
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daily available ground glass and ZrC >2 nanopar- 
tide (NP) doped polyurethane nanocomposites. 
To prepare the nanocomposites we first synthe¬ 
size ZrC >2 NPs using forced hydrolysis [4J] with a 


calcination temperature and time of 600°C and 
1 h respectively. Figure []] shows an SEM im¬ 
age of the NPs at two magnifications and Fig¬ 
ure [2] shows the size distribution of the NPs, 
with the average diameter being 195 ± 32 nm. 
Once synthesized the NPs are then added to an 
equimolar solution of tetraethylene glycol (TEG) 
and poly(hexamethylene diisocyanate) (pHMDI) 
with a 0.1 wt% di-n-butylin dilaurate (DBTDL) 
catalyst, at which point the viscous mixture is 
stirred and centrifuged to remove air bubbles. 
The mixture is finally poured into a circular die 
and left to cure overnight at room temperature. 

The two sample types, ground glass and 
Zr 02 /PU nanocomposite, represent two seperate 
stability and scattering regiemes. The ground 
glass is found to have a scattering length of 970.7 
fjL m and a persistance time on the order of days, 
while the nanocomposite is found to have a scat¬ 
tering length of 11.2 ^urn and a persistance time 
of ~ 90 min. This means that the ground glass 
allows for more stable tests, due to its long per¬ 
sistance time, and larger enhancements due to 
its long scattering length [29|]. On the other 


smaller scattering lengths [291 ]. 

Once the samples are prepared, we use an 
LCOS-SLM based optimal transmission setup to 
compare the different algorithms. The primary 
components of the setup are a Coherent Verdi 
V10 Nd:YV 04 laser, a Boulder Nonlinear Sys¬ 
tems high speed LCOS-SLM, and a Thorlabs 
CMOS camera. See reference [29| for more de¬ 


hand, the nanocomposite allows for tests of the 
algorithms’ noise and decoherence resistance as 
the persistance time is much shorter and noise 
is found to have a larger effect on samples with 


tail. For minimal noise applications we use the 
ground glass sample and operate the Verdi at 
10 W where the laser is most stable (optical 
noise < 0.03% rms) and for measuring the al¬ 
gorithms’ resistance to noise/decoherence we use 
the nanocomposite sample and operate the Verdi 
at 200 rnW where the emission is less stable (op¬ 
tical noise ~ 1%). 

IV. RESULTS AND DISCUSSION 

A. Optimization Speed 

We begin by first measuring the average in¬ 
tensity enhancement for the IA using bins of side 
length b = 16 px, M = 32 phase steps, and a 
ground glass sample. Using the enhancement de¬ 
termined from the IA as the GA stop conditions, 
we measure the number of iterations required for 
each algorithm to reach the set enhancement. 
We use the same bin size for all three algo¬ 
rithms and perform five run averaging of the en¬ 
hancement curves. Figure [3] shows the enhance¬ 
ment as a function of iteration for each algorithm 
with the target enhancement, r/y, marked by the 
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FIG. 1: SEM images of Zr02 nanoparticles at magnifications of 20 OOOx and 200 OOOx. The 
particles are found to be spherical with diameters ranging between 80 nm and 300 nm. 



Diameter (nm) 

FIG. 2: Histogram of NP diameters measured 
from SEM measurements. The mean diameter 
is 195 ± 32 nm. 


the SGA taking 10194 iterations, and the IA tak¬ 
ing 11922 iterations. While the SGA is found to 
be only slightly faster than the IA (1.17x), the 
/iGA is found to be 4.15x faster than the IA, and 
3.55 x faster than the GA for this experimental 
configuration. Further tests with other configu¬ 
rations of samples and bin sizes show that the 
actual speed enhancement varies on the experi¬ 
mental parameters, but for all trials, where the 
three algorithms can reach the same enhance¬ 
ment, the fiGA is the fastest algorithm. This re¬ 
sult is expected given the /uGA’s use of a smaller 
population size 3(1. |31 ]. 


B. Resistance to Noise and Decoherence 


dashed line. From the figure we find that both 
the SGA and /iGA are faster than the iterative 
algorithm, with the //GA taking 2870 iterations, 


In addition to being faster than the IA, the 
GAs are found to be highly resistant to the ef¬ 
fects of both noise and sample decoherence j^ . 
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Number of Iterations (10 ) 


FIG. 3: Enhancement as a function of the 
number of iterations for the IA, SGA, and 
H GA, where the stop condition for the GAs is 
an enhancement of 55. The /zGA is found to be 
the fastest algorithm, with the SGA only being 
slightly faster than the IA. 



FIG. 4: Enhancement as a function of 
normalized iteration using a noisy laser, short 
persistence time sample, and 4096 bins. The IA 
only results in a small enhancement of 3.9, 
while the /iGA obtains an enhancement of 66 
and the SGA an enhancement of 138.5. 


which are c 
IA 


etrimental to the effectiveness of the 


m 


28|, 


29y. To demonstrate the robustness 


of the GAs, we use a Zr02 nanoparticle-doped 


polyuret 
time 


rane sample with a short persistance 


291 ] and operate the probe laser at low 


power where the signal-to-noise ratio is small¬ 
est. First we optimize transmission using the IA 
with a bin length of b = 8 and M = 32 phase 
steps, with the results being poor (shown in Fig¬ 
ure ed with a peak enhancement of 3.9. Next we 
test both the SGA and /i GA with bins of length 
6 = 8 and find drastically improved results from 
the IA with the /zGA reaching an enhancement 


of 66 and the SGA reachin g an enhancement of 
138.5 as shown in Figure [4] 45]. 


Figure 0] also reveals the general functional 
form of the enhancement for each algorithm with 
the IA consisting of many sharp steps, the /zGA 
behaving exponentially, and the SGA behaving 
as a power function. From Figure [4] we find 
that the /zGA rapidly increases and then satu¬ 
rates, with further iterations not adding to the 
enhancement, while the SGA continues to im¬ 
prove until the maximum number of iterations 
is reached. This implies that for the same num¬ 
ber of iterations and bin size the SGA can attain 
















9 


larger enhancements than the /iGA. 


C. /iGA and SGA Comparison 

Given the differences in the performances of 
the /iGA and SGA in the above sections, we mea¬ 
sure the effect of bin size on the two GA’s speed 
and effectiveness. First we measure how bin size 
affects the optimization speed by using the two 
GAs to optimize light transmission through a 
polymer sample to an enhancement of r/ = 43. 
A target enhancement of r/ = 43 is chosen as 
both GAs can achieve that level of enhancement 
easily. Eventually the /iGA is found to saturate, 
while the SGA continues to attain higher en¬ 
hancements (discussed below). This means that 
in order to compare the speed performace of the 
two GA’s we need to use a target enhancment 
that can be reached by both algorithms. Figure[5] 
shows the average number of iterations required 
for the two algorithms to reach the set enhance¬ 
ment. Only bin’s of size b < 16 are shown as 
larger bin sizes are found to be unable to reach 
an enhancement of 43. For all bin sizes tested the 
/iGA is found to be faster than the SGA, with 
a bin size of b = 8 px resulting in the fastest 
optimization for both algorithms. 

In addition to measuring the optimization 
speed of the two GA’s, we also consider the maxi¬ 
mum enhancement achievable given a fixed num¬ 
ber of iterations. Figure [6] shows the peak en¬ 
hancement as a function of the number of bins 
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FIG. 5: Number of iterations required to reach 
an enhancement of q = 43 as a function of the 
total number of bins for the SGA and /iGA. 
Bins with a side length greater than 16 px are 
not shown as they are found to be unable to 
optimize to the target enhancement within a 
reasonable time frame. 


for the GA’s running 12500 iterations. For bins 
of size b > 32, the two GAs are found to pro¬ 
duce the same enhancement within uncertainty. 
However, as the bin size decreases the two algo¬ 
rithms diverge with the SGA producing greater 
enhancements than the /iGA. The divergence 
in the performance of the two GAs occurs as 
the SGA’s larger population is better able to 
explore the solution domain than the /iGA for 
small bin sizes. We also find from Figure [6] that 
for the maximum number of bins (b = 1) the 
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FIG. 6: Peak enhancement as a function of the 
number of SLM bins for the SGA and /zGA 
with a stop condition of 12500 iterations. 
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V. CONCLUSIONS 


We have developed a fi GA for performing op¬ 
timization of light through opaque media. The 
^GA operates based on similar principles to 
previously used genetic algorithms [dj, [ 23 ], but 
differs in the population size, use of elitism, 
crossover technique, and similarity based regen¬ 
eration in place of mutation. These differences 
lead to the /iGA being significantly faster than 
both the SGA and IA. This speed enhancement 
is advantageous for applications in which opti¬ 
mization must occur quickly such as biological 
imaging [liij, authentication 29], and astronom¬ 
ical imaging [hj. For applications where speed 
is less important and large enhancements are de¬ 
sired the SGA is the best option. 


enhancement is much lower than the enhance¬ 
ment for b = 2. This result is slightly counter¬ 
intuitive, as having more bins should allow for 
a more accurate representation of the optimal 
phase mask. However, by increasing the number 
of bins we also increase the solution space size, 
which requires more iterations to effectively ex¬ 
plore. Since we limit the number of iterations, 
we essentially limit the algorithm to searching a 
fraction of the total solution space, which leads 
to smaller enhancements. 


To further enhance the speed of the fiGA we 
are currently working on implementing multi¬ 
threading into the fi GA code in order to take 
advantage of modern computer’s multi-core pro¬ 
cessors. The idea behind using multi-threading 
in the /rGA is to perform overhead calcuations 
(cross over, random number generation, etc.) on 
one thread while a different thread controls the 
SLM and camera. While this speed improvement 
may be small for one iteration, it will compound 
over the use of thousands of iterations to provide 
significant time savings. 
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